Three dimensional super-resolution in metamaterial slab lenses 
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This letter presents a theoretical and experimental study on the viability of obtaining three dimensional super- 
resolution (i.e. resolution overcoming the diffraction limit for all directions in space) by means of metamaterial 
slab lenses. Although the source field cannot be actually reproduced at the back side of the lens with super- 
resolution in all space directions, the matching capabilities of metamaterial slabs does make it possible the 
detection of images with three-dimensional super-resolution. This imaging takes place because of the coupling 
between the evanescent space harmonic components of the field generated at both the source and the detector. 

PACS numbers: 42.30.-d,41.20.Jb,78.70.Gq,78.20.Ci 
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It is well known from the early works of Veselago fl*] that 
a slab made of a left-handed medium will focus the electro- 
magnetic energy coming from a point source to another point 
located at the opposite side of the slab. An experimental con- 
firmation of this focusingof energy has been reported in |2]. 
Subsequent works |3lB13l have shown that, under some cir- 
cumstances, metamaterial lenses can produce images at cer- 
tain planes with a resolution beyond the classical diffraction 
limit, or "super resolution imaging" (SRI). This SRI has been 
attributed to an amplification, inside the lens, of the evanes- 
cent Space Fourier Harmonics (SFHs) coming from the source 
01 ■ In 01 it has been also discussed that this process gives 
rise to fields that decay exponentially from the lens towards 
the image, which causes super-resolution to take place only 
in planes parallel to the slab interfaces. In the direction per- 
pendicular to the lens, a strong decay of the field is observed, 
and thus a three dimensional (3D) picture of the source can- 
not be obtained from the field pattern at the back side of the 
lens. In other words, super-resolution in the transverse direc- 
tions is obtained at the price of a drastic loss of resolution 
in the longitudinal direction. This fact has been corroborated 
by the field measurements recently reported in |4], where the 
field growth from the image plane to the "super lens" can be 
clearly appreciated. Other experimental results also lead to the 
same conclussion, showing that images "of finite depth" Ql 
cannot be directly obtained from field measurements in SRI 
experiments. Although the above facts seem to be well stab- 
lished, new experimental results recently reported by some of 
the authors |7] have suggested the possibility of also obtain- 
ing super-resolution in the longitudinal direction or, in other 
words, three-dimensional super-resolution imaging (3D-SRI). 
In these experiments, a 3D map of a point-like source (i.e. a 
source of sub-wavelength size) was obtained at the image side 
of the lens. The aim of the present letter is then to provide the 
general theory underlying this 3D-SRI. It will be shown that 
3D-SRI of point-like sources is actually a general property of 
metamaterial slabs when the appropriate detection procedure 
is followed. 

The main difference between SRI and more conventional 
imaging processes is that, in the former one, the information 
for the image formation is carried out by evanescent SFHs yfl 



whereas in conventional imaging the information is carried 
out by the propagative SFHs. Since evanescent fields can- 
not carry power, and any image measurement requires some 
power transmission, it is apparent that the dectection of a 
super-resolution image should substantially affect the fields 
around the detector. This intrinsic perturbation of the fields 
in a super-resolution image detection resembles the problem 
of the tunneling effect. As is well known, tunneling of power 
is due to the excitation of a pair of evanescent electromag- 
netic waves, whose interference gives rise to a non-vanishing 
flux of power. "Perfect tunneling" of power in a waveguide 
filled by a metamaterial has been recently reported by some 
of the authors in |8]. In this work it is shown that maximum 
tunneling of power in a setup with identical input and output 
waveguides is achieved when the output is placed at a distance 
from the input equal to that from the source to the image in a 
metamaterial super-lens. This suggests that a similar effect 
could take place in a metamaterial super-lens, provided that a 
detector identical to the source is used for the measurements. 
In the following it will be shown how we can take advantage 
of such effect in order to obtain 3D-SRI with metamaterial 
slabs. 




FIG. 1: Geometry of a metamaterial lens. The source is a loop 
antenna located at z = 0. The detector is also an identical loop an- 
tenna. The lens is formed either by a left-handed metamaterial slab or 
by another planar device that produces SRI through amplification of 
evanescent harmonics. The output antenna can eventually be loaded 
with microwave resistors in order to minimize field perturbation. 
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The present study starts with the canonical problem of the 
formation of images by a left-handed slab of thickness d char- 
acterized by e/eo ~ /^/mo = — 1 + i^, where (5 <C 1 accounts 
for the necessary losses factor to avoid the divergence of field 
integrals | % fllll . In our study, and following an usual 
procedure in microwave SRI experiments, the source will be 
an antenna (specifically, a loop antenna) whose plane is lo- 
cated parallel to the slab interfaces, as is shown in Fig.[n The 
field beyond the lens will be scanned by an output antenna 
that plays the role of detector. The source here employed 
is equivalent to an homogeneous surface distribution of mag- 
netic dipoles inside the loop given by Mg = /qc^"^*, where / 
denotes the amplitude of the imposed time-harmonic current 
in the loop. The computation of the longitudinal magnetic 
field, Hz, at the image side of the slab (z > A + d) can be 
caiTied out by means of the following double inverse Fourier 
transform: 



Hz{x,y,z) = 



G{k,,ky; z)Ms{k,, ky) e^C^^^+fc.?/) , (1) 

where G{kx, ky-, z) and Mg{k^,ky; z) are respectively the 
Fourier transforms of the Green's function of the structure un- 
der study and of the spatial surface distribution of magnetic 
dipoles. After applying the duality principle to expression 
(6) in LI IJ and taking A = d/2, the Fourier-transform of the 
Green's function is found to be 



G(kx 7 ky] Z) 



-kf2/3e- 



-i/3o(2 + d) 



with kf 



t^{(/3/io 



(2) 



kl 



and 13 — y^oj^efi — kf, Po = 
Since Ez = in the present case, the re- 
maining field components can be all deduced from Hz [11]. 

The numerical computation of Q provides the map de- 
picted in Fig.lJlfor the magnitude of Hz in the {x — z) plane 
at the back side of the lens. It should be noted that, given that 
the operation frequency is 3 GHz (Aq = 100 mm) and d and 
A are taken much less than the free-space wavelength, the ray 
model approximation cannot be here employed to obtain the 
fields in the considered region (in other words, we are dealing 
with near fields and therefore in a SRI situation). According 
to previous discussions, Fig.|2]shows that the field magnitude 
has a strong decay along the z direction, so that no informa- 
tion about the location of the source in this direction can be 
extracted from the field pattern. A more detailed picture of 
the field distribution at two planes z = (2 ± d)d is depicted in 
Fig.|3](with 6 set to 0.3), and compared with the field ampli- 
tude near the source at z — ±Sd (both results are identical). 
As is expected from the properties of the metamaterial slab, 
the field distribution at the z = {2 + 6)d plane is almost identi- 
cal to that at the z = Sd plane, thus confirming that the z = Sd 
and the z = {2 + S)d planes are in fact conjugate planes. On 
the contrary, the field distribution at the z = {2 ~ S)d plane 
substantially differs in magnitude from the field distribution 
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FIG. 2: Map of log^Q | Hz \ corresponding to the field generated at the 
back side of the lens by a source loop antenna shown in Fig.Q The 
slab is made of a left-handed medium of thickness d = 4 mm with 
e/eo = m/mo = —1 + iO.OOl, and is separated a distance A = d/2 
from the source. The source is considered, for simplicity, a circular 
loop antenna of radius ri — 5 mm,. The wire radius is of 0.2 mm. 
The operation frequency is 3 GHz. 




FIG. 3: Plot of the logj^Q \Hz \ at different z-planes of Fig.l. The 
slab and the source loop antenna is as in Fig.|2| 



at the z ~ —S plane, near the source. These facts corrobo- 
rate the aforementioned discussions: the super-resolution in 
the transverse directions (the lateral dimensions of the source 
and the image are of about one tenth of the wavelength) is 
compensated by an almost complete loss of resolution in the 
longitudinal direction. 

Next, the problem of the measurement of the image in this 
SRI experiment will be considered. In the microwave range, 
the image detection is performed by measuring the transmis- 
sion coefficient between a source antenna and a receiviiig an- 
tenna, which is scanned in the image side of the lens jHI^I^. 
For simplicity the receiving (or output) antenna will be as- 
sumed to be identical to the input antenna employed as source. 
The transmission coefficient is measured by connecting the 
input antenna to a wave generator via a waveguide, and the 
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FIG. 4: Map of the magnitude in decibels of the transmission coef- 
ficient, 15*21 1, between the input and output antennas of Fig.Q The 
structural parameters are as in Fig.|2| The output antenna is identical 
to the input one. 



output one to a detector through another identical waveguide 
(more details of this measurement setup are reported in 0]). 
In this approach, the metamaterial slab should be viewed as 
a matching device lll2il . whose transmission coefficient, t (or 
in the usual microwave terminology), is given by jjisll 



t = 



.Zu + Zo){Z22 + Zq) 



^12 



(3) 



where Zij are the elements of the impedance matrix for the 
system formed by the two antennas and the left-handed slab, 
and Zo is the characteristic impedance of the input/output 
waveguide. In order to measure a super-resolution image, the 
size of the loop antennas has to be smaller than the free space 
wavelength, which additionally would make the real part of 



Zij (the radiation resistance WlX ) be negligible with respect 
to its imaginary part; namely, Zij ~ —iujLij, where 
is the inductance matrix of the system. The diagonal terms 
of the inductance matrix correspond to the inductances of a 
single loop antenna faced to the left-handed slab. However, 
in the "perfect lens" configuration here considered, the slab 
does not affect the fields around the source and, therefore, 
Lii = L22 = L, where L is the self-inductance of the loops 
in free space. The non-diagonal terms of the inductance ma- 
trix account for the mutual inductance of the loop antennas 
in the presence of the slab, namely L12 = L21 = M. As- 
suming that the waveguides have a low impedance (typically 
this impedance is Zq = 50 17 for microwave measurements), 
it will be found that \Zii\ > Zq. Thus, neglecting second or- 
der terms in Zh/Zq, the transmission coefficient in ^ can be 
approximated as 



2iu;MZo 



(4) 



Note that the only spatial dependence in (0} comes from M 



(since L does not change with the position of the antennas), 
which causes that the transmission coefficient reaches its max- 
imum (\t\ = 1) at those points where M L. Since both the 
input and output antennas are identical, and the field at the 
source plane is reproduced at the image plane (see Fig.|3j, the 
condition M Lis expected to be satisfied around the point 
{x,y,z) = (0,0, 2(i). In consequence, a maximum of the 
transmitted power should be detected at this particular point, 
which thus appears as an "effective focusing point" of the per- 
fect lens. 

In order to show the above expected effects in a practical sit- 
uation, the magnitude of the transmission coefficient has been 
numerically computed for the configuration under study. The 
impedance matrix has been calculated by imposing known 
currents, li = 1,0A, in the loops, and then computing the 
corresponding self and mutual inductances. The computation 
of the self inductance is an standard electromagnetic problem, 
and the mutual inductance is obtained after computing the flux 
of the magnetic field given by across the surface of the 
output antenna. The transmission coefficient is finally deter- 
mined from (|3}, assuming Zq = 5017. In Fig.|4]it is shown 
a map of the computed transmission coefficient for a system 
composed of two identical lossless loop antennas. It can be 
observed that this figure substantially differs from Fig. |2j in 
particular, a clear maximum of \t\ can be observed in Fig.O 
in the neighborhood of the image at (0, 0, 2d). These results 
clearly show the difference between the transmitted power and 
the field distribution in the absence of the output antenna, and 
also how 3D-SRI can be obtained when the appropriate detec- 
tor is used. 

Let us now consider the measurement of the unperturbed 
field. For this purpose, the detector should be designed to 
affect the field distribution as less as possible. It could be 
closely achieved by loading the output antenna with an ad- 
ditional high resistance, which will significantly reduce the 
current induced at the output antenna, as well as its gener- 
ated field. In Fig.|5]the computed values of the transmission 
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FIG. 5: Magnitude of the computed transmission coefficient along 
the Z axis between the input and output antennas of Fig.Qfor differ- 
ent resistances loading the output antenna. The structural parameters 
are as in Fig.|2| 
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coefficient along tlie z-axis of Fig. [T] are shown for different 
resistances, R, loading the output antenna. As is expected, the 
curve for B. = shows a maximum near the location of the 
image, which thus appears as a "focusing point" of the lens, 
at z = 2d. (The small difference between the actual location 
of this maximum and z = 2d can be attributed to the non neg- 
ligible values of the ratio \Zii/Zo\ « 8 in the problem under 
study). As is expected, the curves for the highest values of 
R resemble the behavior of ji/^ | in Fig.|2lfor the unperturbed 
configuration. 

Let us now re-examine the experimental demonstration of 
3D-SRI, recently reported in 1 7] at the light of the above con- 
siderations. Although the underlying physics of the device 
analyzed in |7] is not exactly the same as that of the left- 
handed perfect lens, both systems show the same process of 
amplification of evanescent modes and are equivalent in prac- 
tice for the present purposes. Fig.|6lshows the location of the 
maximum of the transmission coefficient along the z axis of 
the magnetoinductive lens reported in |7], for different values 
of the resistance loading the output loop. The experimental 
setup used to obtain these results is described in detail in Q]. 
The only difference is that, for obtaining the results shown in 
Fig.|6l the output antenna was loaded by different microwave 
resistors. A behavior similar to that of Fig.|5]can be observed, 
thus confirming the proposed theory. 
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FIG. 6: Plots of the measured magnitude of the transmission co- 
efficient between the input and output antennas of Fig. Q when the 
slab is substituted by the magneto inductive lens reported in 01 . 
(Freq=3GHz). 

In summary, SRI in metamaterial super-lenses has been in- 
vestigated. As is well known this imaging is primarily due 
to the amplification of evanescent modes inside the lens. As 
a consequence, super-resolution in the planes parallel to the 
slab is unavoidably compensated by a drastic loss of resolu- 
tion in the direction perpendicular to the slab. However, since 
evanescent modes do not carry power, any physical detection 
(namely, a measurement) of the image at the back side of the 



lens will require the existence of some amount of transmit- 
ted power from the source to the detector, which can signif- 
icantly affect the field distribution. If this last effect is taken 
into account, the secondary fields generated by the presence of 
the detector should be considered. Following this approach, 
it has been shown that the transmission coefficient between 
the source and the detector can be made very high. Provided 
that the appropriate detector is employed (typically a lossless 
output antenna identical to the input antenna), this maximum 
will occur in a neighborhood of the image, thus resulting in 
super-resolution imaging also in the direction perpendicular 
to the slab. However, it should be emphasized that, to obtain 
this 3D-SRI some previous knowledge of the source is neces- 
sary in order to properly design the detector. Thus, a general 
conclusion arises from the analysis: Super-resolution in meta- 
material super-lenses is always incomplete: if the distance 
from the source to the lens is known, the shape and charac- 
teristics of the source can be recovered without uncertainty 
from the analysis of the field at the image side of the lens. 
Conversely, if the shape and characteristics of the source are 
known, it is possible to design an appropriate detector in or- 
der to determine, without uncertainty, the source location in 
three dimensional space from the transmission coefficient be- 
tween the source and the detector. However, it is impossible 
to determine simultaneously, by means of a metamaterial su- 
perlens, both the locations and the shape and characteristics 
of an unknown source wit a resolution overcoming the diffrac- 
tion limit. We feel that the reported analysis and experiments, 
as well as the conclusions arising from them, will be of impor- 
tance in the design of metamaterial super-resolution devices. 
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